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Estimate of average energy of the non-abelian dipole in the instanton medium is calculated. This 
energy for the dipole in colour singlet state escalates linearly with the separation increasing for 
the point-like sources unscreened. And its 'tension' coefficient develops the magnitude pretty 
CD ' similar to one as inferred from the lattice QCD and other model approaches. 

o 

(N . 
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Recent evolution scenario for the QCD matter being created in ultrarelativistic heavy ion collisions 
lO ' proposes interesting explanation of origin of longstanding 'mixed phase' which was found in lattice 
\ calculations of QCD at finite temperatures and, seems, is experimentally observed PQ. A key point to 
^ | account for the phenomenon in Ref. pQ is to justify the presence of colour Coulomb binding interaction 
^ ■ and instanton induced interaction (hence, a lot of bound states available) in rather wide temperature 
q • region above the critical temperature. In spite of many intuitive inputs the ideas of this scenario 
O " (together with some hints at colour singlet states available everywhere in the coupling space in 
Ref. [2]) may help in gaining insight on the QCD dynamics. Our major concern here is to answer 
O question whether the origin of those bound states approaching the critical temperature is instructive 
to understand the nature of ordinary hadrons. 

Since the discovery of (anti-)instantons (the classical solutions to the Yang-Mills equations with 
Ph. nontrivial topological features) the problem of the gauge interaction between instantons is under 
intensive study and the dipole interaction of those pseudo-particles with an external field has already 
been argued in the pioneering papers This very important conclusion was grounded on exam- 
ining the reaction of a single instanton while in a weak slowly varying (constant) external field at 
distances compared to the instanton scale size. We follow that standard approach in considering the 
approximate solution of the Yang-Mills equations as 

a;(x) = a;( x ) + b;( x ) . (i) 

Here the first term denotes the field of a single (anti-)instanton in the singular gauge as we are 
planning to consider an (anti-)instanton ensemble and need to have the non-trivial topology at the 
singularity point 

where p is an arbitrary parameter characterizing the instanton size centered at the coordinate z 
and colour orientation denned by the matrix u>, y = x — z, f\bp,v is the 't Hooft symbol (for anti- 
instanton fj — > 77) and g is the coupling constant of non-abelian field. The second term describes 
an external field. In fact, the result of Ref. [Sj as to the instanton interaction with external field 
(the CDG term) may be obtained in more general way by calculating the corresponding effective 
Lagrangian (without specifying the external field characteristics) For the sake of simplicity we 
limit ourselves here to dealing with SU(2) group only and begin studying an external field originated 
by immovable Euclidean colour point-like source eS 3a and Euclidean colour dipole ±e5 3a which is 
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one of two configurations only providing with the static solutions. Fixing the location of point-like 
sources is gauge invariant unlike specifying their orientation in 'isotopic' space. In order to get rid 
of this 'vice' [S] the solution should be characterized by its energy and total 'isospin' ('isospin' of 
source plus 'isospin' carried by the Yang-Mills field) because these quantities are gauge invariant in 
addition to being conserved. 

With this 'set-up' we could now claim that the results obtained before have maintained the inter- 
acting terms proportional e/g only and in this paper we aim to analyze the contributions proportional 
e 2 which display rather indicative behaviour of asymptotic energy of Euclidean non-abelian point-like 
sources while in the instanton liquid (IL) [Hj and have nothing to do with the additional corrections 
due to the quantum fluctuations about the saddle point. We treat the potentials in their Euclidean 
forms and then the following changes of the field and Euclidean point-like source variables are valid 
B — > iB^, e — > —ie at transition from the Minkowski space. Actually, last variable change is resulted 
from the corresponding transformations of spinor fields ip — > ip, — > 7o — ► 74- It means we are 

in full accordance with electrodynamics where the practical way to have a pithy theory in Euclidean 
space is to make a transition to an imaginary charge. Surely, these Euclidean sources generate the 
fields of the same nature as ones we face taking into account gluon field quantum fluctuations a qu 
around the classical instanton solutions A = A inst + a qu . 

For Euclidean point-like source and dipole in colourless state we take the Coulomb solutions of 
the classical Yang-Mills equations in the following forms 



B a (x) = (0,5 a3 V),<p = ^- 




(3) 
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respectively, with z e being a coordinate of point-like source in peace and z±, z 2 as coordinates 
characterizing a dipole. As known, the field originated by particle source in 4d-space develops the 
cone-like shape with the edge being situated just at the particle creation. Then getting away from 
that point the field becomes fully developed and steady in the area neighboring to it. In distant area 
where the field penetrates into the vacuum it should be described by the retarded solution obeying, 
in particular, the Lorentz gauge d^B^ = (which is also valid for Eq.(j2J)). Here we are interested in 
studying the interactions in the neighboring field area and the cylinder-symmetrical Coulomb field 
(which is £4 independent as in Eq.Q) might be its relevant image in 4d-space. 
Thus, the field strength tensor is given as 

G% = dpAZ - d v A« + g e ahc A\A c v , (4) 

where e abc is the entirely asymmetric tensor and for the field superposition of Eq. (0) it can be given 
by 

G% = G%(B) + GIM) + G%(A, B) , 

G%(A,B) = ge abc (B b ^A: + AlB:) , (5) 
if G^^A) and G a ^ v {B) are defined as in Eq.Q. Then gluon field strength tensor squared reads 

g% g; u = g; u (b) g; u (b) + g; u (a) g%(A) + g%(a, b) g; u (a, b) + 

+ 2 G^{B) G%(A) + 2 G* U (B) G* U (A, B) + 2 G* U (A) G%(A, B) . 6 

The various terms of Eq.© provide the contributions of different kinds to the total action of full 
initial system of point-like sources and pseudo-particle 

S = Jdx(± G%, G%, + 3 ;Al) = S e (B) +(3 + S mt . (7) 
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So, the first and second terms of Eq.© account for the source self-energy S e (for the dipole it should 
be supplemented by the Coulomb potential of interacting sources) and the single instanton action 
(3 = 87r 2 /g 2 . At small distances the first term is regularized by introducing the source 'size' as in 
classical electrodynamics. These terms are proportional to e 2 and g~ 2 , respectively and, clearly, out 
of our interest here. The fourth and last terms of Eq.© together with the term j^A® of Eq.(|7|l 
account for the contribution to the interacting part of the action 5^ proportional to e/g. And 
the third (repulsive) term of Eq.© only leads to the contribution proportional to e 2 as the fifth 
term is equal zero owing to the gauge choice. Collecting those non-zero contributions we have after 
performing the calculations the following result 



S m t = - Vk4i ^ k h + (£) 2 J + (£) K kl luV . (8) 

The exact form of I{ (where the CDG term is absorbed) is unnecessary in this paper because going 
to relate an external field to the stochastic field of fluctuations (the IL model) we have to average 
over the colour orientation of (anti-)instanton. It leads to disappearance of the dipole contribution 
and the first nonvanishing correction, as known, comes from a condensate term (but in next order). 
Two other terms result from the 'mixed' component of field strength tensor which looks like 
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G%(A, B) = 2^ ^ fj kia p WT ^ { 



-,2 



GUA, B)=l4- E a3C U Ck fj kia % -, J- T s I , — r-r - t— 

4jV ; 4?r y (y + p ) \\y + Ai| \y 



(9) 



for the point-like source and for the field of colourless dipole, respectively, and where A = z — z e , 
Ai i2 = z — Zip- The other contributions to G°L, (A, B) are absent because of the gauge used. In 
order to handle the further formulae easily it is practical to introduce new dimensionless coordinates 
as x/p — > x. Then for the single source the function J and the tensor K take the following forms 



J =2 dy 



2y 2 -y 
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Kh — 2 J dy 



y A {y 2 + if\y + ^\ 2 ' 
y k yi 1 

"7" W+WW+Af ' 



and can not be integrated in the elementary functions. Fortunately, we need to know their asymptotic 
values at A — > oo only 

5tt 2 1 

/ ~ 

2 A 2 ' 

and for the components of the 2-nd rank tensor 

Kij = 5 i:j K x + AjAj K 2 , 

we have 



K x ~ —-2 , K 2 ~ . 



2 A 2 

Clearly, the 'mixed' component of field strength tensor is of purely non-abelian origin but its con- 
tribution to the action of whole system (point-like sources and pseudo-particle) takes the form of 
self-interacting Euclidean source ~ e 2 although it is descended from the instanton field and field 
generated by sources. It seems, this simple but still amazing fact was not explored properly. 

Now we are trying to analyze the pseudo-particle behaviour in the field of Euclidean non-abelian 
source and develop the perturbative description related to the pseudo-particle itself. Realizing such 
a program we 'compel' the pseudo-particle parameters to be the functions of 'outside influence', i.e. 
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putting p — > R(x,z), uj — > Q(x,z). These new fields-parameters are calculated within the multipole 
expansion and then, for example, for the (anti-)instanton size we have 



R in {x,z) = p + c^y^ + c^ y fl y v + ... , \y\ < D 

y y vu (10) 

R out {x,z) = p + d M -| + d^/-f-2+..., \y\> D , 

v y y 



Similar expressions could be written down for the (anti-)instanton orientation in the colour space 
Q(x, z) with D being a certain parameter fixing the radius of sphere where the multipole expansion 
growing with the distance increasing should be changed for the decreasing one being a result of 
deformation regularity constraint imposed. Then the coefficients of multipole expansion c M , c^, . . . 
and dp, d^ u , . . . are the functions of external impact. It turns out this approach allows us to trace 
the evolution of approximate solution for the deformed (crumpled) (anti-)instanton as a function of 
distance and, moreover, to suggest the self-consistent description of pseudo-particle and point-like 
source fields within the approximate solution of Eq.Q (our paper has been completed recently). The 
role of deformation fields was found to be essential in the local vicinity of source (at distances of 
order < 2p) and results in considerable growth of the pseudoparticle action. 

Proceeding to the calculation of average energy of point-like source immersed into IL we have to 
remember that one should work with the characteristic configuration which is the superposition of 
instanton and anti-instanton fields also supplemented by the source field B®(x), i.e. 
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Al(x) = Bl{x)+Y, A l&n) » (11) 



i=l 



where ji = {pi, Z{, ioi) are the parameters describing the i-th (anti-)instanton. This would-be su- 
perposition ansatz is a solution to the Yang-Mills equations only in the limit of infinite separa- 
tion but approximates well enough the dominant configuration saturating the functional integral 
[3]. The IL density at large distances from source practically coincides with its asymptotic value 
n(A) ~ n e /3 ~' s ~ n Q because the action of any pseudo-particle here is approximately equal (3. The 
quantity we are interested in should be defined by averaging S over the pseudo-particle positions and 
their colour orientations (taking all pseudo-particles of the same size at the moment) as 

dzi 



_1 (12) 
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e 2 1 „ „ „ „ fdA/e\ 2 /, K , 
An a 

here L is a formal upper integration limit, V = L 3 X^ defines the IL volume, X4 is an upper bound of 
the 'time' integration, N denotes the total number of pseudo-particles and a is a source 'size' value 
(on strong interaction scale, of course). Taking the asymptotic functions of J and K and returning 
to the dimensional variables for the moment one instanton contribution to the average action can be 
written down in the following form 

P 2 1 /V N fi7T 3 P 2 

(S) - f - X 4 + £ /3 L 3 X 4 + ^- e - f L X 4 , (13) 
47r a V V p g 

where p is the mean size of (anti-)instanton in IL. The result is given in the form with the common 
factor X4 extracted because our concern now is the (anti-)instanton behaviour in the source field 
background where the field is fully developed and the solution possesses the scaling property at any 
time slice. In the limit N, V — > 00 and at the IL density n = N/V fixed this result becomes 

E^^- + nPL* + n%- 2 p 2 L. (14) 
4vra /3 g 2 
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The last term here looks like a correction to the gluon condensate (the second term). However, 
this contribution linearly increasing with L and proportional to e 2 has different physical meaning of 
additional contribution to the source self-energy 

67T 3 6 2 

E ~a L , a = n — — 5 P 2 ■ ( 15 ) 
P 9 

From the view point of the CDG 'ideology' [3 j an appearance of such a term is quite understandable 
if we remember that for long-range Coulomb solutions Eq.(j3J) the tensor G a falls like r~ 2 at large 
distances (similar to merons and regular instantons), markedly slower than one for the instanton 
(~ r~ 4 ). Besides, the colour 'magnetic field' behaviour for our configuration (of time slice) with 
point-like source looks as 

Hi = -5° 



ra3 X k 



\X 
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bringing the result to the formal analogy with the Polyakov's model of vacuum at a fixed time albeit 
its physical meaning here is much more obscure and rather far from the interpretation discussed in 

The terms next to leading order could be obtained by a cluster decomposition of the exponential 
with the interacting term derived, i.e. Eq.(|5|) 

(exp(-S M )) UM = ex P ( E > ( 16 ) 

where the indecies uz just imply aforementioned averaging over the pseudo-particle positions and 
their colour orientations and (Si) = ((Si)), (S1S2) = (S1XS2) + ((S1S2)), . . . Then, for example, the 
second order contribution generated by dipole interaction Jj looks like 



(„-!)£. 

and results in (due to nontrivial dielectic IL properties [3]) the regular correction to the average IL 
energy (fTlj) as 

N e 2 4vr 3 

IT V" 3, 

Now dealing with the field of colour dipole in the colour singlet state we are able to demonstrate 
the resembling reaction of IL once more. The contribution to average energy E is defined at large 
distances by the integral very similar to that for the configuration with one single source 

dA, ( X - 2 , A I A , + — 7 ) ^ (1> 



A^ (/) = -^T f- (17) 




Ax |Ai||A 2 | 

When the source separation I = \z\ — z 2 \ is going to zero, the field disappears and the final result 
should be zero. There are only two parameters to operate with the integral, they are L and I. The 
dimensional analysis reveals that the integral is the linear function of both but the /-dependence only 
obeys the requirement of integral petering out at I — > 0. It is easy to receive the equation 

+ L = I 

for determining the coefficient (the contributions of three integrals are shown separately here). Finally 
the average dipole energy for the colour singlet (c.s.) state reads 




E ~ a I , (c.s.) , (19) 
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The 'tension' value a for the IL characteristic parameters ~ *t where R is the mean separation 

R 3 

of pseudo-particles, n = i?~ 4 , (3 ~ 12, p ~ 1 GeV^ 1 [S],[Z! comes about a ~ 0.6 GeV/fm (if one 
takes for the brief estimate for the source intensity e ~ g) . This result looks fully relevant in view of 
the qualitative character of estimates which IL is able to provide with. Moreover, such a value is in 
reasonable agreement with the estimates extracted from the potential models for heavy quarkonia, 
for example. If one intends to explore the magnitudes like (Si nt e~ Smt ) (which could model an 
effect of suppressing pseudo-particle contribution in the source vicinity) in numerical calculations 
it becomes evident the linearly increasing behaviour starts to form at A/p ~ 3 — 4. Thus, the 
asymptotic estimate obtained shows the Euclidean source energy in IL gives the major contribution 
to the generating functional in quasi-classical approximation if all the coupling constant are frozen 
at the scale of mean instanton size p. 

The linearity of static potential obtained delivers, as we believe one interesting message. The 
energy growth with the distance increasing hints strongly at a 'great difficulty' of bringing the 
Euclidean colour source into IL because the source mass (the additional contribution received should 
be treated just in this way) is unboundedly increasing if the Debye screening effects do not enter the 
play. It can be shown also that corresponding integral in the color triplet (c.t.) state leads to the 
result 

with color triplet mean energy 

E~a(AL-l), (c.t.), (21) 

what clearly demonstrates again it is difficult to have the states with open colour in IL. Besides, 
it encourages that studying the instanton correlations and collective degrees of freedom induced by 
them could be the promising way to grasp the full QCD dynamics. 

If one strives to go away from the approximation of specifying the source orientation in the 
'isotopic' space then the natural generalization of the results obtained takes the form ('sources of 
large intensity') 

^ = (PxPi)L + 2 (P 1 P 2 ) (hP2)L = {P x + P 2 f L - {PiP 2 ) I , 

and eP\, eP 2 are considered as an intensity of the sources with P = (P 1 , P 2 , P 3 ) denoting the vector 
of the unit normalization \P\ = (p a p a y/ 2 = \ m colour space (the details of two body problem 
analysis for arbitrary source orientations in colour space could be found in [E],[II]). Surely, the result 
demonstrates again that it is very difficult to reveal the states with open colour in IL and besides 
that there is another small parameter in the problem. Indeed, the deviations of vector-sources from 
the anti-parallel orientations may not be large. 

Another point which is worth of a discussion here concerns the Wilson loop behaviour. The 
absence of the area law for that in [Sj is explained by rapid decay (~ A~ 4 ) of the corresponding 
correlators in IL and inadequate contribution of large size instantons. However, what we are able 
to estimate now emphasizes an essential impact of the part of average G 2 originated by the inter- 
ference term of colour source field and stochastic (instanton) component when one determines such 
observables as the interaction energy of the probe (infinitely massive) particles. This quantity devel- 
oping the linear increase at large distances transforms into the constant gluon condensate while in 
the stochastic background only. (Even in vacuum the energy required to separate two probe colour 
charges from each other is finite if there are the light quarks.) Perhaps, the physical picture could 
be made more transparent if other observables are considered. 

The estimate of both configurations contribution to the generating functional (see Ref.jH] for 
notation) 

Z = T,Tn ft / WPi) Cn c e~ s , (22) 

N iV - i=l J 
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can be performed with the variational maximum principle. Then the mean energy of Euclidean 
sources appears in the exponential factor and may be interpreted as the suppression factor for the 
states with open colour 

Z > e~ EXi , 

here we omitted the condensate contribution together with residual Coulomb component. Calculating 
the Wilson loop with this weight we find the area law behaviour for the colour dipole in the 'isosinglet' 
state but have no transpaparent physical picture to take it as a signal of forming the string solution. 

We believe the result above shows strong interaction vacuum picture based on the localized lumps 
of topological charges (which is highly successful in describing many nonperturbative QCD phenom- 
ena and hadron phenomenology but at the same time was failling in treating a confinement) still have 
enough credibility to become a cogent theory of QCD vacuum providing all-round understanding of 
strong interaction dynamics on the same footing in lieu of a set of individual mechanisms. 
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